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After simplifying and improving the non- Gaussian formalism we developed in previous work, we de- 
rive a quantitative expression for the three-point correlator (bispectrum) of the curvature perturba- 
tion in general multiple-field inflation models. Our result describes the evolution of non-Gaussianity 
on superhorizon scales caused by the nonlinear influence of isocurvature perturbations on the adia- 
batic perturbation during inflation. We then study a simple quadratic two-field potential and find 
that when slow roll breaks down and the field trajectory changes direction in field space, the non- 
Gaussianity can become large. However, for the simple models studied to date, the magnitude of 
this non-Gaussianity decays away after the isocurvature mode is converted into the adiabatic mode. 



INTRODUCTION 

The comparison of inflationary predictions with in- 
creasingly precise observations has become a major fo- 
cus of cosmology over the past decade. Observations 
of the cosmic microwave background (CMB), notably 
by WMAP and the forthcoming Planck satellite, pro- 
vide accurate measurements of the CMB temperature 
anisotropies. If this experimental effort is to place any 
further constraints on inflation, it should be matched by a 
corresponding effort to improve the accuracy of theoreti- 
cal predictions. To date, the study of linearized perturba- 
tions around a homogeneous background has provided an 
excellent first approximation. However, gravity is inher- 
ently nonlinear so it is important to investigate whether 
higher-order effects are relevant for future high-precision 
observations. 

Any nonlincarity will manifest itself as non- 
Gaussianity in the correlations or other measures 
of the stochastic properties of the CMB anisotropies, 
a prime measure being the three-point correlator or 
bispectrum. The general expectation has been that 
non-Gaussianity generated during inflation will be very 
small and unobservable. In this paper, summarizing 
and extending results from our two earlier papers 0, 0| , 
we quantitatively study this question for multifield 
inflation and apply our results to models with two scalar 
fields, evaluating our expressions numerically. We find 
that non-Gaussianity evolves due to the influence of 
isocurvature perturbations on the adiabatic perturbation 
and can become large when slow roll breaks down and 
the field trajectory changes direction in field space 
- this happens when the heavier field rolls down to 
reach its minimum value and the lighter field takes 
over. However, in the simple cases studied to date, this 
large non-Gaussianity is temporary and decays after 
the turn of the field trajectory, as the conversion of the 
isocurvature modes is completed. 

The effort to characterize non-Gaussianity from infla- 
tion has attracted a great deal of interest because it may 



open up a new observational window on the early uni- 
verse (see e.g. the review 0] and more recent papers 
0, HI). This effort has been successful in investigating 
non-Gaussianity in a variety of specific regimes, but in 
this paper we focus on the important influence of the 
isocurvature modes on the adiabatic mode during infla- 
tion. First, we introduce the basic formalism using our 
key gradient variable Q to characterise nonlinear inhomo- 
geneities in the long- wavelength approximation. Next, 
wc identify a constraint between the first isocurvature 
mode and the adiabatic time derivative, which allows us 
to reduce the 2n nonlinear evolution equations to an el- 
egant 2n— 1 system (with n fields). In the two-field case, 
we then perturbatively expand these equations to sec- 
ond order and we obtain an integral expression for the 
three-point correlator in terms of linear solutions and 
background slow-roll parameters (at no stage requiring 
these to be small). A numerical investigation of a two- 
field quadratic model then demonstrates how a poten- 
tially large bispectrum emerges and decays. 

NONLINEAR INHOMOGENEOUS QUANTITIES 

We begin by summarising and extending the formal- 
ism developed in [2j to compute the nonlinear evolution 
on superhorizon scales during multifield inflation. Space- 
time is described by the long-wavelength metric 

ds 2 = -N 2 (t, x)dt 2 + a 2 (t, x)dx 2 , (I) 

with the lapse function N(t,x), defining the time slic- 
ing, and the scale factor a(t, x) both dependent on space 
and time. The key variable describing the full nonlinear 
inhomogeneities is Q 

C = S ml di lna - (e mA d^ A ) , (2) 
V 2e 

with k 2 = 8irG = 8"7r/mpj and 4> A the scalar fields of 
the multifield inflation model; e and e A n are defined be- 
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low. Note that we allow for the possibility of a nontriv- 
ial field metric Gab, used to raise and lower indices in 
field space. The Einstein summation convention is as- 
sumed implicitly throughout this paper. The nonlinear 
quantity Q n transforms as a scalar under changes of time 
slicing on long wavelengths Q ; under the transformation 
X^ = (t, x) -> X^ = (T(t, x),x), 



C(t,x)=C(T,x). 



(3) 



Here, the spatial part of the transformation is indepen- 
dent of time and can be taken to be trivial [2j. In the 
single-field case and when linearised, ^ is just the spa- 
tial gradient of the well-known curvature perturbation £ 
in the literature. Note that our variable has been 
further studied in a covariant framework in Q and its 
relationship to perturbation theory elaborated. 

We have chosen an orthonormal set of basis vectors 
in field space, which are defined from the field ve- 
locity, e A = II A /II, and successively higher-order time 
derivatives via an orthogonalization process 0, @] • Here 
II A = <ft A /N is the proper field velocity, with length 77. 
For example, e A = {8^-e A e 1B )V t B B /| . . . |, withX> t the 
time derivative covariant with respect to the metric Gab 
in field space. A similar orthogonalization procedure for 
the special case of two fields and a flat field metric was 
independently proposed in Q . This basis has the advan- 
tages of immediately separating the adiabatic component 
Q from the isocurvature modes C z m {jn > 2) and making 
the origin of non-Gaussianity more transparent. We also 
define the local slow-roll parameters 0, V 



e(t,x) = 



k 2 B 2 
277 2 
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r,W(t,x) = -3- 
+ 3(e-77 ll )<5 m i-3 ? 7 J 



n A v A 
bb 2 ' 



(4) 
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where V(<p A ) is the potential of the model, with its co- 
variant field derivatives Va = 9aV and V-ab = T>bOaV. 
Furthermore, V mn = e^e^V-AB and 77 ee a/{Na) is the 
Hubble parameter. We have chosen the number of c- 
folds t = In a (i.e. TV 77 = 1) as the time variable, which 
simplifies the equations. We emphasize that we are not 
making any slow-roll approximations, the slow-roll pa- 
rameters are just short-hand notation and can be large. 

In the definitions ([J}, e and 7/H are the effectively single- 
field slow-roll parameters. The truly multiple-field slow- 
roll parameter r/ 1 - is proportional to the component of 
the field acceleration perpendicular to the field velocity; 
i.e. it is non-zero only when the trajectory in field space 
makes a turn (at least for a flat field manifold; if the 
manifold is curved even a straight trajectory can have 
77 7^ due to the connection terms). To physically in- 
terpret x we rewrite it as 377 2 x = V22 — Vn + 6B 2 e, valid 
to leading order in slow roll. Thus, \ is related to the dif- 
ference between the effective mass perpendicular to the 



field trajectory and the one parallel to it, with a cor- 
rection proportional to 77 2 e, which is actually dominant 
initially. In this language SH 2 ^ — —V21 (leading-order 
slow roll) is proportional to the off-diagonal mass V21. 
The £ m are second-order slow-roll parameters. 

In the field and Einstein equations were reformu- 
lated as an equation of motion for the inhomogeneous 
variable Q n . We combine Q a and its time derivative 
0™ ee dtiCi 1 ) m t° a single vector, 



w ia = (Ci ) &i 7 Ci , @i , Ci ) #i 



(5) 



with T denoting the transpose. (The reason for the 
change of notation with respect to |l[ will become clear 
below.) Then the nonlinear evolution equation, valid on 
superhorizon scales, becomes 



( (t, x) + B a p(t, x)wip(t, x) = 0, 



(6) 



where the index a runs over twice the number of fields. 
This equation is exact within the long-wavelength ap- 
proximation where second-order spatial gradients are ne- 
glected and is valid beyond slow roll. The full expression 
for the matrix B a p can be found in for two fields on 
a flat manifold in the iV77 = 1 gauge it is exactly 
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e + 27/H -67/- 1 - Aer/ 1 - - 2£ J 
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-2t/ j 
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where 



3x + 2e 2 +467/H +£H 3 + e + 2r]«J 

(7) 



Since the matrix B is spa- 
(j> A (t, x), and n A (t,x) [see 
has to be closed with rela- 



£1 and ^ ee £ 2 . 
tially dependent via H(t, x) 
(J4j)] , the nonlinear system ([£ 
tions expressing these quantities in terms of Wi a . Again, 
the exact constraint equations, derived from the Einstein 
equations, can be found in but in the two-field case 
they become 



8,, In 77 = e(t, e mA d^ A = -— C\ (8) 



et-DiBA 



V2e 

K 

K \ 

It turns out that there is an exact algebraic relation be- 
tween 9} and C 2 , which allows us to reduce the order of 
the system by one. The derivation of this relation goes as 
follows. Taking the m = 1 component of the T>iB A con- 



straint and multiplying by — n/(H\'2e), we find that the 
left-hand side equals —di In 77. Taking the time deriva- 
tive of this equation and comparing with the equation of 
motion for B\ [the a = 2 component of ([6])] , using the fact 
that d t In 77 = 7/" as can be seen from the field equation, 
V t n A = -377 A - G AB V. B /B, we obtain the announced 
exact long- wavelength result: 



JLy»2 



2r/ x C. 



(9) 
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A linear version of this result can be found in [8| , but here 
we show it is valid nonlinearly as well. It is the general- 
ization of the well-known result [13] that C is constant in 
single-field inflation, showing the influence of the isocur- 
vature (e2) mode on the adiabatic (ei) mode in multiple- 
field inflation. Note the result is valid for an arbitrary 
number of fields, not just for two fields. 

Hence we can now define a vector vi a without the 6\ 
component: 

v ia = (d,c?,olc!,el...) T , (io) 

satisfying our key nonlinear long-wavelength evolution 
equation [2| with the constraint imposed, 

v ia (t, x) + A ab (t, x)v lb {t, x) = 0, (11) 

where the index a runs over 2x (number of fields) — 1. For 
two fields on a flat manifold in the NH = 1 gauge the 
matrix A is given by 

/0 -2r) x \ 

A = -11. 

\0 3x + 2e 2 + 4er/" +4(77^)2 3 + e + 2r?H/ 

(12) 

For a curved manifold with a nontrivial field metric the 
term — {2e / K 2 )e2 Rabc d&i ef should be added to the 



A = 



with 




where V Xmn = (V2~e/ K )efege%V A Bc/(3H 2 ). Note that 
in the NH = 1 gauge we are using here, the exact A a b c 
only contains terms up to next-to-leading (third) order 
in slow roll, no higher-order terms exist. 

The source term on the right-hand side of the first- 
order equation (| 13() is not present in our original long- 
wavelength equation (llip and requires further explana- 
tion. In this form equation (I13p is nothing more than the 
full linear perturbation equation for Q n rewritten after 
smoothing over short wavelengths with a window func- 



A 32 component, with R A BCD the curvature tensor of the 
field manifold. Again, we stress that no slow-roll approx- 
imation has been made. 



PERTURB ATIVE EXPANSION 

To solve the master equation (JTTJ) analytically, we ex- 
pand the system as an infinite hierarchy of linear pertur- 
bation equations with known source terms at each order. 
To second order we obtain 

** (i . , +^s ) w«a ) - as) 

C + ^k? = (*,*)«}?, d4) 

where v la = v\ 1 } + vf}, and 

A ab (t,x) = A%i+A$ = AW+d-WfrA^W 

= ^(t)+A^(t)v^(t,x). (15) 

Here we have denoted = d~ 2 d l v^ and diA a b is com- 
puted using ([8]) and ([4]) . For two fields on a flat manifold 
the exact result for A a b c is 



tion W(kR) with smoothing length R = c/(aH) and c a 
constant of the order of a few. 1 The window function is 



In the gauge t = ln(aH ) used in p], Q] , R is a function of time 
only (which was the reason we chose that gauge for a stochas- 
tic numerical implementation). In the present gauge, t = In a, 
which is more convenient for super-horizon calculations, this is 
no longer true. However, the spatial dependence of H is of higher 
order in slow roll [see |(H}], so that it can be neglected if slow roll 
holds during horizon crossing. 



2e^ -4t?H?7- l + 2£- l 
-6 X ~ 2677 11 -2(7 7 H) 2 - 2(r/ ± ) 2 
-6 - 






-2e 2 - 467/H +2(r/H) 2 - 2{ V ± ) 2 -2^' 
-4er/- L - 2£- L 
-2r/- L 



(16) 



- 6(t/ ± ) 2 - 3e X - 4e 3 - 10e 2 7/H - 2e(r/H) 2 - 6e(7/ ± ) 2 + 8r/H (r/^) 2 - 3ef " - 6?/^ + |(7 m - V 22 i) s 
-12er/- L - 6t/V + 1277- L x ™ 6eV + 4(t/II)V + 4(r/- L ) 3 - 4e£- L - 2r/H£- L + |(F 2U - V222) 

6?/- 1 - - 2er/- L + 47/H77- 1 - - 2£- L 

I 
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chosen such that: (a) W(kR) -> for k > 1/R and (b) 
W(kR) — > 1 for fc <C 1/-R. The first condition ensures 
that short wavelengths are cut out, making our long- 
wavelength approximations applicable, while the second 
ensures that at sufficiently late times the solution of (fT3"|) 
does not depend on the exact shape of W, simply being 
the appropriate linear solution on long wavelengths. The 
source bi a can then be expressed in terms of the linear 
mode function solutions x2rn{k,t) as 

&& = / ^^W(fc)xW(fc)at n (fc)iA :i e ifc -+c.c.. (18) 

Because of the conditions (a) and (b) described above, 
the time derivative W(fci?) peaks around a time just af- 
ter horizon crossing and the term 6^ can be seen as 
providing the initial conditions for the long-wavelength 
linear evolution of each mode k. Given this behaviour of 
W(kR), the matrix X am in (fT8|) . containing the Fourier 
solutions from linear theory for £ and 9, must be known 
around horizon crossing. These linear solutions can be 
found directly numerically, but a general analytic solu- 
tion is known in the slow-roll approximation [8J]. In the 
two-field case, the slow-roll horizon-crossing solution is 
simply [H 

X am = " , ^= \ 1 | • (19) 

The quantum creation (a} n ) and conjugate annihilation 
operators in (|18p satisfy the usual commutation rela- 
tions. 2 

In the second-order equation (Tl4"|) , the source term on 
the right-hand side represents nonlinear corrections due 
to long-wavelength evolution. One might be concerned 
that any nonlinearity produced before and during hori- 
zon crossing (b\J) is ignored. However, recent work ex- 
aminin g b oth tree- level [12j (see also [l|) and quantum 
effects [13[ indicates that any such nonlinearity will be 
suppressed by small slow-roll factors and lead to non- 
Gaussianity of the order of that produced by single-field 
inflation. Hence any significantly larger non-Gaussianity 
will have to come from the subsequent nonlinear evolu- 
tion on superhorizon scales (which is a purely multiple- 
field effect), which is why we concentrate on this regime 
in our treatment. 

Equations (fT3")) and (TT4")) (as well as any higher order 
equations), together with the initial condition Vi a (t — > 



2 in [H 01 we used a stochastic description of the source term, nec- 
essary for a complete numerical implementation. Here we can re- 
tain the standard quantum description of [Tl[. This leads to dif- 
ferent \/2 factors in some intermediate steps compared with [lj . 



— oo) = implied by the formalism, can be solved using 
a single Green's function G ab (t,t'), satisfying 

— G ab (t,t')+A^(t)G cb (t,t') = 0, KmG ab (t,t') = S ab . 

The solution is the time integral of G ab contracted with 
the terms on the right-hand side of (|13| and (fT4|) . Hence, 
the superhorizon linear mode functions v2m are given by 

«$t(M) = / dt' G ab (t,t')W(k,t')X^(k,t'), (21) 

J — OO 

with 

v$(t,x) = I ^^t,W(fc,i)at n (fc)ifc 4 e ifc - + c.c.. 

(22) 

On superhorizon scales, the solution (f2"Tj) is equivalent 
to that from the full linear perturbation equation, so we 
could have obtained it from a direct numerical calcula- 
tion rather than projecting forward our analytic solution 
at horizon crossing as we did here. In particular, the 
integral in (|2ip is independent of the exact form of the 
window function W for times sufficiently long after hori- 
zon crossing (in practice this means just a few e- folds). 
We note that a convenient choice for explicit calculations 
is to employ a simple step function ('top hat') as win- 
dow function, so that VV(fc,i) = 6(kR/^/2 — 1), with 
R = (c/ H) exp(— i) the smoothing length, which is dif- 
ferent from the Gaussian window function used in [l[. 
Finally, then, at second order we obtain the expression 

= - / dt'G ab (t,t')A bcd (t')v^(t',x)v^(t',x), 

(23) 

which is again independent of the exact choice for W. 

BISPECTRUM 

So far we have used time slices for which NH = 1, be- 
cause it simplifies superhorizon calculations, that is, time 
slices in which the expansion of the universe is homoge- 
neous (t — In a) . On such time slices 

d = --^(e 1A d^ A ) = ^d iP , (24) 

where p is the energy density. [To derive this, one uses 
©, Q, and H 2 = n 2 p/3.} One can show that Q is 
a total gradient in the single-field case. However, for 
multiple fields this is not the case and hence Q cannot 
be written as the gradient of some scalar. To connect 
with scalar observables it is more convenient to go to 
uniform energy density time slices where dip = [6(. On 
such time slices we simply have 

(I =d i \na = d i a, (25) 



5 



Denoting these new time slices by T — t + At and ex- 
panding the right-hand side of ([3]) to second order we 
have (suppressing the x dependence) 



Ci (2)1 (*) = (l 2)1 (t)-Atd t C 



(i)i 



(26) 



So no correction appears at linear order, which is as it 
should be, since Q is gauge-invariant at linear order. At 
second order there is a correction proportional to the 
time derivative of Q, which is zero in the single-field case. 
The quantity At is the time difference between uniform 
expansion and uniform energy density time slices, and is 
evaluated by comparing a in the two gauges: 



a(t) = a(T,x) = a(t + At,x) 
= a(t, x) + At d t a(t, x) + 



(27) 



which, expanded to first order, becomes 
a (0) (t) = a (0) (t) + a™ (t, x) + At^ (t, x)d t a^ (t). (28) 



Hence we find that = and 



At 



ftaC°) 



-c 



(1)1 



(29) 



(NH = 1 => dtcS^ = 1) so that we can rewrite (|26[) as 

9,a( 2 )=cf )1 = Cf 1 +2r ? x C (1)1 C| 1)2 . (30) 

Using our formulae, we can now write down the power 
spectrum and bispectrum for a — d~ 2 d l Q. For the 
power spectrum we have 

(aaf\k,t)=vV {k , t )v{V(k,t), (31) 

with xj^l in ([2lTl assumed to be real [as in 



(|19[)]. By combining the different permutations of 
(a^ (xi)a^ (x2)a^ (x 3 )) containing the linear and 
second-order adiabatic solutions, we arrive at our main 
result for the connected bispectrum: 



(aady (fci, k 2 , k 3 ,t) 



(32) 



2tt) 3 5 3 (£ s k s ) f(h,k 2 ) + f(h, k 3 ) + /(fc 2 , k 3 ) 



with 



f(k, k') = L mn (k, k' , t)v^ [k, t)v§ (k',t)+k~ k' (33) 



and 



L mn (k, k', t) = r,Ht)4n(k, t)v$ (k', t) 



(34) 



d*'GUM04°l(OC(M0^(fc',*') 



(i)/ 



This expression, valid if the non-Gaussianity produced at 
horizon crossing can be neglected, gives the bispectrum in 



f_NL (eq . mom . } 
f_NL (squeezed) 
-f]par 
?7perp 



/ / 



TA \ 



FIG. 1: The non-Gaussianity parameter /nl (|36l) for the two 
momentum configurations described in the text together with 
the slow-roll parameters — r)" and r] J ', plotted as a function 
of time in the model (|35p with initial conditions <f>\ = (j>2 = 
13k -1 and mass ratio m2/m\ — 12. 



terms of background quantities and linear perturbation 
quantities at horizon crossing. It is exact: no slow-roll 
approximation has been used. Since Q = dia is a total 
gradient, there are no nonlocal terms involving inverse 
powers of momenta. 3 In a forthcoming paper [l4| we will 
further work out this expression for the bispectrum an- 
alytically. Here we finish with a numerical investigation 
of a specific model. 



NUMERICAL INVESTIGATION OF A 
TWO-FIELD QUADRATIC POTENTIAL 

We investigate a simple model with two fields and a 
quadratic potential: 



(35) 



choosing m\ = 10 k (the overall mass magnitude can 
be freely adjusted to fix the amplitude of the power spec- 
trum) and initial conditions 0i = 02 = 13k" 1 for a total 
of about 85 e-folds of inflation. The higher the mass ratio 
m 2 jm\ is chosen, the sharper a turning of the trajectory 
in field space this model has during the last 60 e-folds, 
which is reflected by a peak in T) , possibly followed by 
additional oscillations. 

In figure [1] we show the result for the non-Gaussianity 



3 ln Eli where we did not make the switch from 1 to the total 

gradient Q , the result contained factors of the form k ■ (k + 
k')/\k+ k'\ 2 . One can easily relate that result to this one by 
noting that k- (fc + k') = k'\ 2 + \{k 2 - k' 2 ), where the 

second term disappears after the correction from switching to 
dip = time slices has been included. 
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parameter /nl, defined as [l[ 

(aaa)(k 1} k2,k 3 ) 



NL 



[(aa)ki (aa)k 2 + (fe <-> fe) + (&i <-> fe)] /3 



(36) 



[where the (27r) 3 and the ^-function are removed from 
the expression for the bispectrum] , in the case of a mass 
ratio mzjrnx = 12. It is plotted as a function of time 
during the few e-folds of the corner-turning in field space, 
together with rj» and rj- 1 . We show the results for two 
momentum configurations: (a) all three momenta equal, 
with horizon crossing 60 e-folds before the end of infla- 
tion, and (b) two momenta equal, crossing the horizon 60 
e-folds before the end of inflation, with the third momen- 
tum crossing the horizon 20 e-folds earlier (squeezed con- 
figuration). We see that a large /nl is produced initially 
but, in this model, is completely erased by the time r] 1 - 
has its first zero. 4 The dependence on different momen- 
tum configurations can be inferred from (fT§|) and (f2"Tj) . 
The main momentum dependence is the overall factor of 
fc~ 3 / 2 of vim! however, if other momentum dependence 
were absent, these factors would cancel exactly in the 
definition of /nl- The remaining dependence on scale 
comes from the slow-roll suppressed dependence of X am 
on time and the fact that VV(fc, t) peaks at different times 
for different modes. Because of this slow-roll suppression 
the momentum dependence is rather weak; figure [1] shows 
a variation of less than a factor 2 between two very dif- 
ferent configurations. 

The large non-Gaussianity of the type described above 
might be accessible in models where inflation ends dur- 
ing the corner-turning, such as hybrid inflation, or in 
which residual isocurvature modes persist. We have also 
investigated quartic potentials, and find that large, but 
temporary, non-Gaussianity emerges in a similar fashion. 
It may be that persistent large non-Gaussianity is easier 
to achieve in less symmetric potentials. Analytic calcu- 
lations in 14| , as well as further numerical investigations 
that are currently underway, should settle these issues. 



DISCUSSION 

In this paper we extended and improved the non- 
Gaussian formalism developed in our earlier work. We 
derived and made use of a long-wavelength constraint 
relation (O, which simplifies the system of equations de- 
rived in []], Q]. We also clarified the issue of different 
time slicings at second order in order to properly con- 
nect to scalar observables. Using this formalism we have 



4 Because of subtle cancellations between terms of different order 
in slow roll, the leading-order slow-roll estimate in Q (for the 
case of mass ratio 9) overestimates the non-Gaussianity substan- 
tially, although it does provide a fair estimate for the maximum 
value reached during inflation. 



provided a general, exact expression for the bispectrum 
in multifield inflation due to superhorizon effects, equa- 
tions (|52"j) - (|M|) , which provides a powerful way to com- 
pute non-Gaussianity for an arbitrary inflationary model 
up to the end of inflation. The computation is technically 
straightforward since it only requires the background so- 
lution and the linear perturbation solution at horizon 
crossing to be known, either analytically or numerically. 

A numerical investigation of the simplest two-field 
quadratic example showed that non-Gaussianity can be- 
come large, as parametrized by a value of /nl of order 
unity, whereas single-field inflation generically predicts 
a result two orders of magnitude smaller [15|. In this 
model a temporary breakdown of slow roll is required 
somewhere during the last 40-50 e-folds of inflation (i.e. 
some time after the observable scales have left the hori- 
zon so that the scalar spectral index will not be too far 
from unity). This happens for mass ratios larger than 
about 10 when the trajectory of the scalar fields driv- 
ing inflation turns a corner, as parametrized by r\ . The 
large non-Gaussianity is produced by the nonlinear influ- 
ence of the isocurvature mode on the adiabatic mode but 
it also decays away as the isocurvature mode disappears. 
The exact behaviour of /nl as a function of time dur- 
ing inflation is complicated by subtle cancellations and a 
full explanation is left for a future publication [3. We 
find that there is only a weak dependence on the specific 
momentum configuration. 

We note that earlier work in [l6| emphasised the im- 
portance of the relation between non-Gaussianity and 
isocurvature modes, although their models relied on non- 
linear interacting potentials to generate higher-order cor- 
relators. Related work in [B| with the l SN' formalism, 
using the same quadratic two-field potential, produced 
results which are qualitatively similar, that is, with large 
non-Gaussianity emerging temporarily. Numerically the 
SN approach evolves neighbouring inflationary trajecto- 
ries through a multidimensional field space, with appro- 
priate derivatives estimated from these. In contrast, our 
approach relies more directly on integrals over the linear 
perturbation and background solutions which, in prin- 
ciple, should be more accurate. A direct quantitative 
comparison and an analysis of the relative merits of the 
two approaches deserves further investigation. 

Attempts to incorporate inflation within the frame- 
work of a more fundamental physical theory like string 
theory or extensions of the Standard Model invariably 
lead to inflationary models with many scalar fields. Our 
findings indicate that non-Gaussianity can emerge nat- 
urally in these models, although its perdurance seems 
closely linked to that of the isocurvature modes. It is 
of great interest to investigate scenarios in which this 
non-Gaussianity might survive to late times because its 
imprint in the CMB may be observable in forthcoming 
experiments. 
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